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Abstract. Autocorrelation functions (ACFs) are studied for a sample of 16 long gamma-ray bursts (GRBs) with known 
redshift z, that were observed by the BATSE and Konus experiments. When corrected for cosmic time dilation, the ACF shows 
a bimodal distribution. A narrow width class (11 bursts) has at half-maximum a mean width t'^ = 1.6 s with a relative dispersion 
of ~ 32%, while a broad width class (5 bursts) has r'^ = 7.5 s with a ~ 4% dispersion. The separation between the two mean 
values is highly significant (> 7cr). This temporal property could be used on the large existing database of GRBs with unknown 
redshift. The broad width set shows a very good linear correlation between width at half-maximum and (l+z), with a correlation 
coefficient R = 0.995 and a probability of chance alignment < 0.0004. The potential application of this correlation to cosmology 
studies is discussed, using it in combination with recently proposed luminosity indicators. 
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1. Introduction 

The knowledge of time scales and source distances are es- 
sential for the physical understanding of astronomical phe- 
nome na. From the first detections in 1969 by Vela satel- 
htes dKlebesadel. Strong. & 01sonlll973l) . until the launch of 
BeppoSAX in 1997, the distance scale of gamma-ray bursts 
(GRBs) remained unsettled. This mission provided arc-minute 
localization, leading to the discovery of a fading emission to- 
wards lower energy bands, the so-called afterglows. Thereafter, 
burst redshifts z have been determined from spectroscopic anal- 
ysis of the afterglows or, in some cases, of their associated host 
galaxies, proving that at least long-duration bursts are at cos- 
mological distances. So far, the redshift of no short-du ration 
burst has been clearly determined (although see lKulkarni et alJ 
l2002h . In this paper only the class of long GRBs will be con- 
sidered (i.e., those with time duration > 2 s). 

Up to date, more than 30 burst redshifts have been spec- 
troscopically measured thanks to immediate follow-up obser- 
vations. On the other hand, there is a wealth of data from 
thousands of GRBs for which the redshift is unknown. Most 
of these were detected by the Burst and Transient Source 
Experiment (BATSE). Other important motivations to find 
a redshift estimator based only on the gamma-ray prompt 
emission are the lack of optical counterparts in some cases 
(the so-called dark afterglows), and the difficulty of spec- 
troscopically determine redshifts beyond z - 5 due to the 
Lyman alpha absorption. In recent years, two empirical re- 
lations have been discovered to estimate the luminosity dis- 
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tance exclusively from the analysis of the gamma emission. 
One relates the isotropic luminosity wit h the ti me lag between 
different energy channels (Norris et al.i l200cil) , and the other 
with a variability parameter of the light curve jReichart et alJ 
l200lli . Both luminosity correlations can be used to derive 
luminosity distances and, assuming some specific cosmol- 
ogy, the coiTesponding redshifts. Thus, from these correla- 
tions it has been possible to estimate GRB lun iinosity func- 
tions and demographic distributions (see, e. g., iNorrisI E'002t 
iLlovd-Ronning. Frver. & Ramirez-RuiJ2002l) . These first esti- 
mations indicate that the GRB population may peak at redshift 
z ~ 10, being then ideal probes of the early universe. However, 
the luminosity functions derived in these works predict source 
counts A^(> P), as a function of photon flux P, that differ sig- 
nificantly from the observed one (Schmidt 2003). Much better 
calibration of these empirical relations is needed, and that will 
only be possible with a much larger number of independent 
redshift determinations covering a broader z range. 

Individual power density spectra (PDS) of GRB are in gen- 
eral very diverse, but the longest bursts show power-law spec- 
tra extended over two frequency decades. Shorter bursts also 
display this property by averaging the PDSs o f a large sample 
( iBeloborodov. Stern. & S vensson 19981 l2000l) . This underly- 
ing power-law behavior indicates the absence of any preferred 
time scale. The autocorrelation function (ACF) is the Fourier 
transform of the PDS, therefore it contains in principle the same 
information that can be visualized in a different way. The ACF 
gives a measure of the correlation between different points in 
the light curve that are separated by a given time lag. Various 
efforts have been made using these data analysis tools to find 
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a temporal characteristic that might correlate with the redshift 
and, e.g., Chang, Yoon, & Choi (2002) have found a weak cor- 
rela tion be t ween the power-law index of the PDS and z- See 
also lAttei'al (l2003l) for a proposed redshift indicator. 

In this paper it will be shown that the ACF can be used 
to define characteristic times that strongly correlate with the 
redshift. In § |2 the data selection and the use of the ACF are 
described. Next in §|3lit is shown that the ACF corrected for 
time dilation effects has a bimodal distribution, and that this 
property could be used to construct an empirical relation to es- 
timate z. Finally, the results and their possible applications are 
discussed in §|3 



2. Data and methods 

This work is mainly based on data taken by BATS E on board 
the C ompton Gamma-Ray Observatory (CGRO; Fishman et alJ 
It consisted of eight modules placed on each corner of 
the satellite, giving full sky coverage. Each module had two 
types of detectors: the Large Area Detector (LAD) and the 
Spectroscopy Detector (SD). The former had a larger collecting 
area and from it the CGRO Science Support Center (GROSSC) 
provided the so-called concatenated 64 ms burst data, which 
is a concatenation of the three standard BATSE data types 
DISCLA, PREB, and DISCSC. All three data types have four 
energy channels (approximately 25-55, 55-1 10, 1 10-320, and 
> 320 keV). The DISCLA data is a continuous stream of 1 .024 
s and the PREB data covers the 2.048 s prior to the trigger time 
at 64 ms resolution, both types obtained from the 8 LADs. They 
have been scaled to overlap the DISCSC 64 ms burst data, that 
was gathered by the triggered LADs (usually the four closer to 
the line of sight). This combined data format was used when 
available, since the concatenated pre-burst data allows a better 
estimation of the background. 

All BATSE bursts with known redshift z were considered 
for study ' . In some cases, like GRB 980326 and GRB 980613, 
the data are incomplete or were not recorded. For burst GRB 
970828 the DISCSC data are incomplete; but it was possible 
to derive data with the same characteristics from the MER data 
type, binning up the 16 energy channel into 4 DISCSC-like 
energy channels. For GRB 000131 the given DISCSC data are 
unevenly sampled and it had to be uniformly binned into 1 .024 
s time resolution. The BATSE sample total 1 1 cases. 

To improve statistics, it was also considered GRB data that 
are p ublicly available fro m other experiments. The Konus mis- 
sion ("Apteka r et alJl99 '?) publishes GRB light curves of 64 ms 
resolution within an energy band of 50-200 keV. At the time 
of this publication, there were 25 Konus bursts with known 
redshift. But the collecting area of this experiment is about 
20 times smaller than the one on BATSE and consequently, 
in most cases, the signal is too weak for the purposes of this 
analysis. A total of 5 bursts were selected for this study (read 
further discussion below). 

The INTEGRAL mission (launched in October 2002) 
makes public all count time histories of the bursts detected 



by the anti-coincidence shield of its gamma-ray spectrometer 
(SPI-ACS). They have a time resolution of 50 ms and a non- 
sharp lower energy threshold at about 80 keV dKienlin et alJ 
COOL) . So far, the only detected burst with known z is GRB 
030329, and it was also detected by Konus. Therefore, these 
data were used here mainly for comparative purposes. 

The au tocorrelation function of GRBs was fi rst studied by 
iLink et al ( 1993) and l ater on by, e.g.jFenimor e et all d 19951) 
and lBeloborodov et'aP 1I2OOOI) . Following their notation, from 
a uniformly sampled count history with AT" time resolution and 
time bins, let m, be the total observed counts at bin /. Also 
let bi be the corresponding background level and c, = m,- - bi 
the net counts. The discrete ACF as a function of the time lag 
T = kAT is 



Air) = 2 
1=0 



l,...,N- 1 



(1) 



and A(0) - 1 for k - 0. Here the periodic boundary condi- 
tions (c, - ci+n) are assumed. The normalization constant A,, is 
defined as 



N-l 



(2) 



1=0 



' See J. Greiner web page for a compilation of all GRBs with known 
redshift at _http://www.mpe.mpg.de/~jcg/grbgen.htnil^ 



The normalization makes the ACF of each burst fluence in- 
dependent. The term m, in Eq. |2l subtracts the contribution of 
the uncorrected noise assuming that it follows Poisson statis- 
tics. This is not the case for the pre-trigger data of the con- 
catenated BATSE data type and they were excluded after the 
background fitting. Their contribution to the ACF is negligible 
but they would affect the estimation of A„. This should also 
be considered when doing power density spectral analysis. For 
practical reasons, the actual calculation of Eq.^was done us- 
ing a Fast Fourier Transform (EFT) routine, i.e., squaring the 
absolute value of the transformed gives the power density spec- 
trum, and taking the inverse transform of it gives the ACF. Zero 
padding of the time series was used to avoid the artifacts pro- 
duced by the periodic boundary condition. 

The background estimation was done by fitting with up to 
second order polynomial the pre- and post-burst data, that was 
judged by visual inspection to be inactive. This is particularly 
critical for weak bursts. Unfortunately, the Konus GRB light 
curves that are publicly available have a fixed duration of 100 
s, with no pre-burst data, and sometimes not even post-burst 
data. Only a few cases are sufficiently bright and have long 
post-burst data to allow a reliable estimation of the ACF. The 
problem was studied using numerical simulations and it be- 
came clear that for most of the Konus set the systematic errors 
introduced by the background estimation are the main source of 
uncertainty. Figure^shows comparisons of the ACFs of bursts 
for which there are data from both the Konus and BATSE ex- 
periments. As reported by Fen imore et a l. (1995), the ACF of 
GRBs narrows at higher energies. The best match was obtained 
using the sum of the BATSE energy channels 2 and 3, cover- 
ing a similar energy range that the corresponding Konus data. 
Note that the agreement will depend mainly on having a simi- 
lar lower-end energy limit, since there are more counts at lower 
energies and the ACF is a quadratic function of the number 
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Fig. 1. Comparison of the ACFs of 6 GRBs obtained using data from two different experiments. Solid lines: Konus 64 ms data in 
the 50-200 ke V energy band; gray lines: BATSE 64 ms data in the 55-320 keV energy range. There is sufficiently good agreement 
for the bright bursts, when the resuhs are not very sensitive to the background estimation. GRB 971214 is considerably dimmer 
than the others (see the text for discussion). 



of counts. The Konus weak case GRB 971214 illustrates how 
a poor estimation of the background affects the ACF calcula- 
tion. On the other hand in the strong case GRB 990123, even 
with a short post-burst data tail to fit the background, the differ- 
ence between the ACFs is acceptable for the purposes of this 
work. Guided by this comparison, the selection criteria for the 
Konus cases were set, requiring peak count rates larger than 
3000 counts s"' and post-burst data. These criteria are met by 
all bursts shown in Fig.^except the first, and by 5 other cases 
not observed by BATSE that were then added to the sample. 
Among these last cases is the bright GRB 030329 that was also 
observed by INTEGRAL, and Fig.|2lshows the good agreement 
between the ACFs derived using the two different instruments 
data. Tabled summarizes in its four first columns the adopted 
sample of GRBs, the instrument source, the estimated redshift 
z, and the corresponding reference. 

3. Results 

In Fig. the ACFs of all the GRB sample are shown. One 
can see that at different heights the width of the ACFs has a 
fairly uniform distribution, with the half-width ranging, e.g., 
at half-maximum between 2.5-20 s. Figure |2' shows the lo- 
cal autocorrelation function A(t'), where the cosmic time di- 
lation has been removed, and r' = t/(1 + z) is the corrected 
time lag. The width of the different ACFs shows now, par- 
ticularly around the half-maximum level, a bimodal distribu- 
tion with a clear gap between two sets. A broad width set of 
5 bursts and a narrow width set of 1 1 bursts. The redshift of 
GRB 980329 is only known to be in the range z = 2.0-3.9 
( iLamb. Castander. & Reicha rt"l999). Thus, for Fig.Et an av- 
erage value z - 3 was chosen, but in any case for the given 
z range its ACF will lie within the other narrow width bursts. 
This burst was used here to show the bimodality but will be 
excluded from the following calculations. For the 1 1 BATSE 
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Fig. 2. Two ACF functions of GRB 030329. Solid line: Konus 
64 ms data in the 50-200 keV energy band; gray line: 
INTEGRAL 50 ms data with a soft low energy cut-off at ~ 80 
keV. The high secondary peak of the ACF reflects the two 
bright pulse structure of the light curve. 



bursts the local ACF distribution was analyzed at different en- 
ergy channels. Although narrower at larger energies, the ACF 
shows the same clear bimodal distribution in all channels. 

To study the distribution of the local ACF, e.g., to esti- 
mate statistical moments at different lags, one cannot simply 
add the time series since now all of them have time bins of 
different sizes. To overcome this problem, the logarithm of 
each discrete ACF was approximated by a polynomial function 
/(t) ^ In A(t). a high degree polynomial was used (typically 
~ 12) to match the data within the range of the random fluc- 
tuations up to time lags of 10 and 30 s for narrow and broad 
cases respectively. These ranges were chosen to well cover the 
central peaks of the ACFs down to the 0. 1 level. Using these 
functions, the mean and the sample standard deviation s were 
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Fig. 3. a) Autocorrelation functions of the 16 bursts sample. BATSE and Konus cases are shown in gray and solid lines respec- 
tively, b) Local ACFs, where time dilation due to cosmic effect has been corrected, being t' = t/(1 H- z). GRB 980329 is shown 
with dashed lines assuming a redshift z = 3. 



calculated for the two sets. Since the sample size n is small in 
both cases, the standard deviation cr was estimated as 



using a correction for low number statistics assuming normal 
distributions. Figurel^shows the mean curves and the ±lcr re- 
gion around them. At half-maximum A{T'g) = 0.5, the half- 
width is = (1.6 + 0.5) s and t'^ - (7.5 ± 0.3) s for narrow 
and broad width sets respectively. Hence, at half-maximum the 
distributions do not overlap even at the 7cr level. The proba- 
bility of such separation into two sets by chance, having an 
underlying unimodal distribution, was estimated numerically. 
The overall distribution was assumed uniform in a given range 
(most favorable case). Considering the same total number of 
points (15), and asking conservatively the size of the smallest 
set to be > 4, the probability p of obtaining two sets with a dif- 
ference between their means > 7cr is p < 6 x 10"^. Therefore, 
the alternative hypothesis of having a bimodal distribution can 
be accepted at a highly significant confidence level. 

The small range of the broad width distribution around the 
half-maximum level is particularly interesting, because it rep- 
resents a relative dispersion of 4%, while the relative dispersion 
of the narrow width is 32%. This means that if we had a way 
to know a priori at what width class a burst belongs, we would 
be able to estimate with the same corresponding uncertainties 
the time dilation factor I + z, and therefore the redshift in prac- 
tice only when z > 1 . The width t„ was calculated fitting the 
logarithm of the ACF in the range 0.4 < A(t) < 0.6 with a sec- 
ond degree polynomial. In the last two columns of Tabled the 
obtained values for and t'^ - To/(1 + z) are listed. Figure |5] 
shows To versus I + z for both width classes. As expected, the 
broad width set shows a very good correlation, with a linear 
correlation coefficient R = 0.995 and a probability of chance 
alignment p < 0.0004. The corresponding values for the nar- 
row width set are R = 0.809 and p < 0.005 respectively. Notice 
that GRB980425, which has been associated with SN 1998bw, 
belongs to the broad width set. This burst was considered an 
outlier in the studies of the lag and variability luminosity cor- 
relations when modeling the data with a single power-law, al- 
though its inclusion supports the general trend in both cases. 
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Fig. 4. Mean value of the local ACF for narrow and broad width 
bursts (solid lines). The Icr region about the mean is also shown 
(gray lines). At half-maximum A(t') - 0.5, t' - (1.6 + 0.5) s 
and r' = (7.5 + 0.3) s for narrow and broad widths respectively. 



4. Discussion and conclusions 

The average PDS of bursts shows an overall power-law be- 
havior, indicative of a self-similar underlying process where 
there are no preferred timescales. If this is the case, then the 
width of the ACF is related with the low-frequency cut-off 
of the PDS, which is du e to the finite duration of the burst 
("Beloborodov et al.'2000). As it was mentioned in §n in prin- 
ciple the information given by the ACF and the PDS is the 
same. In practice, since they express this information differ- 
ently, they are affected by noise and window effects in different 
ways. It would be difficult to make a good estimation of the 
low-frequency cut-off in the PDS due to the large statistical 
fluctuations. On the other hand, the width at half-maximum of 
the ACF gives a robust measure. 

In this analysis the ACFs of bursts were only corrected for 
the cosmic time dilation. However, since the detectors are sen- 
sitive over a finite energy band, effects due to the shift in en- 
ergy s hould also be present. Studing a set of 45 bright long 
bursts, iFenimore et al.l ( Il995h found that the full-width W of 
the average ACF (at the e"-^ level) depends on the energy E 
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Table 1. Sample of 16 GRBs with known redshift. The 6 
columns give the name of the GRB, the instrument, the mea- 
sured redshift z, the corresponding reference, the ACF half- 
width at half-maximum Tq, and the width corrected for time 
dilation t' . 



GRB 


Instrument 


z 


Ref. T,(s) 
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2 


15.33 
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3 


8.02 


1 81 




R AT<^P /Tt'oniic 
ocj/ js.unus 


'X -t- 1 


4 


5.96 


1 AQ 


980425 


BATSE 


0.0085 


5 


7.62 


7.56 


980703 


BATSE 


0.966 


6 


14.15 


7.19 


990123 


BATSE/Konus 


1.600 


7 


19.81 


7.62 


990506 


BATSE/Konus 


1.3066 


8 


3.83 


1.66 


990510 


BATSE/Konus 


1.619 


9 


2.54 


0.97 


991208 


Konus 


0.7055 


10 


3.67 


2.15 


991216 


BATSE/Konus 


1.02 


11 


3.80 


1.88 


000131 


BATSE 


4.500 


12 


5.77 


1.05 


000210 


Konus 


0.846 


13 


2.43 


1.31 


010222 


Konus 


1.477 


14 


3.68 


1.48 


011121 


Konus 


0.362 


15 


9.98 


7.32 


030329 


Konus/INTEGRAL 0.1685 


16 


2.56 


2.19 



9 
10 
11 
12 
13 
14 
15 
16 



Metzgeret^jl997h. 



Dioreovski et alJ<200ll) . 

Kulkamietal.'('l998t) . 

Lamb, Castander, & ReichartI jl999h . 
Tinnev et a l. 1 1998). 



Dioreovski et aLU99'i) . 
KulkarmetalJh999|) . 



Bloom et al. "7200 




Beuermann et alJ ^1999 
Dodonov et al. (1999, 
Vreeswij^^L 1.19991) . 
Andersen et aljj200(i) . 



Piro et al 
Jha et al, 



7200: 



200D. 



Gamavich et alj<2q03h . 
GreineretalJ<2003h . 



as WiE) cc This narrowing of the ACF should partially 

counteract the time stretching since for large redshifts the en- 
ergy window of the instrument will see photons emitted at 
higher mean energies. Furthermore, due to the trigger threshold 
bursts detected at high redshifts are more luminous. There are 

indications that the pulse width, and therefore the ACF width, 

1 ' 

correlates with the luminosity (Lee, Bloo m. & Petros ian 2000). 

One should consider also that because of the energy shift bursts 
at high redshifts are detected at earlier stages. If the local ACF 
has an approximately constant width t', (for each width class) 
these effects should produce a deviation from linearity in Fig.|5] 
Since no important deviation is observed, the net combined ef- 
fect must be small. To explore how sensitive our results are to 
such effects we will assume now that the local width of the 
ACF is given by t'^ = To/(1 + z)'^", where the index a takes 
into account additional redshift dependencies. Figure |5] shows 
the relative dispersion of the width for each set as a function 



20 



15 



t9 10 



broad width t'„ a 
narrow width t'. ■ 



1 2 3 4 5 6 

1+Z 

Fig. 5. Correlation between the width at half-maximum To and 
the dilation factor 1 + z. Bursts are classified into narrow and 
broad width r'^ cases, as seen by a local observer, and they are 
marked by squares and triangles respectively. GRB 980329 is 
indicated with a hollow square; its redshift is only known to 
be within the shown z range. A linear fit of the form f(x) = ax 
gives fl„ - (1.42±0.15)s andfli - (7.54±0. 10) s for narrow and 
broad sets. The slopes are approximately equal to the average 
for each set. Note that since the slope ratio is 5.3, for a given 
To each correlation predicts very different 1 -i- z values. 



0.5 

0.4 

^ 0.3 
"... 

e 

0.2 
Q.1 



— bnoai^ 

\ narrow 

- s / 

" ^ - 
















-0.4 -0.2 





0.2 


0.4 



8 

'1 



Fig. 6. Variation of the dispersion of the ACF width with addi- 
tional redshift dependencies. The local values are calculated as- 
suming = To/(l H-z)'^", where a is an unknown correction in- 
dex. The left ordinate shows the relative dispersion of the broad 
(dotted line) and narrow (dashed line) sets. Dispersion minima 
occur at fl <s; 1 in both cases. The right ordinate shows for 
the solid line the difference between the mean values of each 
set (jUfo, fin) in number of total standard deviations (cTb + cr„). 
The gap between sets remains above 3cr over a large range 
(-0.4 < a < 0.3), indicating the robustness of the bimodal- 
ity result to additional corrections. 



of a. The dispersion minima occur at small a index values in 



both cases, with a 



(h) 



-0.05 ± 0.05 and a 



(«) 



-0.13 ±0.23 



for the broad and narrow width sets respectively. The difference 
between the mean values of each set versus a is also shown in 
Fig-El It peaks at a = -0.05 with 8.8cr, where now cr is the total 
standard deviation. Noteworthy, the gap between sets remains 
larger than 3cr over a large range (-0.4 < a < 0.3), indicating 
how robust the bimodality result is to any additional correction. 
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The practical use of the proposed empirical relations re- 
quires a criterion to decide to what width class a burst belongs. 
For a width Tg < 7.5 s the two classes do not overlap (see 
Fig-E}. On the other hand, if extrapolations are valid, t„ » 20 
s implies for the narrow width class unrealistically large z- No 
burst spectral characteristic was found to correlate with the 
width class. The same was true for the luminosity and total en- 
ergy release, but again larger samples should be studied to be 
conclusive. The mentioned luminosity correlations could give 
a first z estimation to determine the width class, and then it 
will be possible using the ACF to obtain a second better and 
independent estimation. Since the two classes are separated by 
more than a factor 5, the selection should not depend in prac- 
tice on the assumed cosmological parameters, and therefore 
the ACF width-redshift correlation could be used in addition 
to constrain them. 

To investigate the combined use of these correla- 
tions, the ACFs of a B ATSE burst sample studied by 
iFenimore & Ramirez-RuizI (|2000.) were calculated. They esti- 
mated the redshifts of 220 GRBs using a power-law fit of the 
luminosity L versus the variability V based on only 7 bursts. 
Therefore they obtained the best-fit parameters for thek model 
L oc V" with large uncertainties. Furthermore, the best-fit pa- 
rameters vary considerably if GRB 980425 is included in the 
calibration, and this burst was excluded in the derivation of the 
published redshifts. Figure shows the logarithmic distribu- 
tion of the observed ACF widths Tq derived for this sample 
(dashed-line) and the distribution of the widths corrected 
for time dilation {solid-line). In Fig. those cases where the 
iterative method to determine z diverged have been excluded, 
as well as bursts with observed ACF widths t,, < 2.4 s, but 
our conclusions do not depend on this selection. While the un- 
corrected distribution seems unimodal, the corrected one ap- 
pears bimodal, in qualitative agreement with the results of §|3] 
The probability p that a statistical fluctuation could produce 
this feature was estimated. Assuming an underlying log-normal 
distribution, a conservative estimation gives p < 0.02. The 
distribution is considerably broad, but that is to be expected 
given the large spread found in the luminosity-variability cor- 
relation. The median values for the two width subsets are ~0.7 
s and ~4.0 s, both approximately a factor 2 smaller than the 
mean values of the distributions shown in Fig. |3 However, 
this discrepancy can be accounted considering the uncertain- 
ties. Based on the an alysis of 20 bursts with known redshifts 
[Reichart et al.' (2001) estimated a — 3.3^,) ^. In particular, the 
exclusion of GRB 980425 will overestimate the exponent a and 
consequently the redshifts, given smaller r'g corrected widths. 

The ACF width-redshift correlations described inhere 
will need to be confirmed by a larger statistical sample. 
Additionally, the lag and variability luminosity correlations 
need to be known for a larger redshift range to avoid uncertain 
extrapolations. The close agreement between the ACFs using 
data from the past mission BATSE and the presently operating 
Konus and INTEGRAL will allowed to continue improving the 
statistic of this work. In combination with the luminosity cor- 
relations we should be able to construct a GRB-based Hubble 
diagram (i.e., a luminosity distance versus z plot) for high z, fol- 
lowing a procedure similar to that of iSchaeferA2003.) . To obtain 
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Fig. 7. Distribution of the ACF width at half-maximum for a 
sample of 170 BATSE GRBs. Dashed-line: logarithmic his- 
togram of the width To (in seconds). The distribution seems 
unimodal. Solid-line: The width corrected for cosmic time dila- 
tion t'o, using estimated redshifts derived from the luminosity- 
variability correlation ( Fenimore & Rarnirez-Ruiz 2000). The 
distribution appears now bimodal, with median values approx- 
imately 0.7 s and 4.0 s for each subset. 

such diagram would have important implications in cosmology 
studies. Ongoing efforts in this direction will be presented in 
the near future. 
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